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I describe harmonic-oscillator-based effective theory (HO BET) and explore the 
extent to which the effects of excluded higher-energy oscillator shells can be 
represented by a contact-gradient expansion in next-to-next-to-leading order 
(NNLO). I find the expansion can be very successful provided the energy de- 
pendence of the effective interaction, connected with missing long-wavelength 
physics associated with low-energy breakup channels, is taken into account. I 
discuss a modification that removes operator mixing from HOBET, simplifying 
the task of determining the parameters of an NNLO interaction. 
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1. Introduction 

Often the problem of calculating long-wavelength nuclear observables - 
binding energies, radii, or responses to low-momentum probes - is formu- 
lated in terms of pointlike, nonrelativistic nucleons interacting through a po- 
tential. To solve this problem theorists have developed both nuclear models, 
which are not systematically improvable, and exact numerical techniques, 
such as fermion Monte Carlo. Because the nuclear many-body problem is so 
difficult - one must simultaneously deal with anomalously large NN scat- 
tering lengths and a potential that has a short-range, strongly repulsive 
core - exact approaches are numerically challenging, so far limited to the 
lighter nuclei within the Is and lower Ip shells. The Argonne theory group 
has been one of the main developers of such exact methods.^ 

However, effective theory (ET) offers an alternative, a method that lim- 
its the numerical difficulty of a calculation by restricting it to a finite Hilbert 
space (the P- or "included" -space), while correcting the bare Hamiltonian 
H (and other operators) for the effects of the Q- or "excluded" -space. Calcu- 
lations using the effective Hamiltonian H'^^^ within P reproduce the results 
using H within P + Q, over the domain of overlap. That is, the effects of Q 
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on P-space calculations are absorbed into P[H'^^^ — H)P. There may exist 
some systematic expansion - perhaps having to do with the shorter range 
of interactions in Q - that simplifies the determination P{H^f^ — H)P?'^ 

One interesting challenge for ET is the case of a P-space basis of har- 
monic oscillator (HO) Slater determinants. This is a special basis for nuclear 
physics because of center-of-mass separability: if all Slater determinants 
containing up to n oscillator quanta are retained, H*^^^ will be translation- 
ally invariant (assuming H is). Such bases are also important because of 
powerful shell-model (SM) techniques that have been developed for inter- 
ative diagonalization and for evaluating inclusive responses. The larger P 
can be made, the smaller the effects of W^^^ — H. 

There arc two common approaches to the ET problem. One is the de- 
termination of P[H^-ff — H)P from a given H known throughout P + Q, 
a problem that appears naively to be no less difficult than the original 
P + Q diagonalization of H. However, this may not be the case if H is 
somehow simpler when acting in Q. For example, if Q contains primarily 
high- momentum (short-distance) interactions, then H'^^ — H might have 
a cluster expansion: it becomes increasingly unlikely to have m nucleons in 
close proximity, as m increases (e.g., a maximum of four nucleons can be 
in a relative s-state) . Thus one could approximate the full scattering series 
in Q by successive two-body, three-body, etc., terms, with the expectation 
that this scries will converge quickly with increasing nucleon number. This 
would explain why simple two-nucleon ladder sums - g-matrices - have 
been somewhat successful as effective interactions'* (however, see Ref. 5,6). 

The second approach is that usually taken in effective field theories,^ 
determining H'^^^ phcnomcnologically. This is the "eliminate the middle- 
man" approach: H itself is an effective interaction, parameterized in order 
to reproduce NN scattering and other data up to some energy. So why 
go to the extra work of this intermediate stage between QCD and SM-like 
spaces? This alternative approach begins with PHP, the long-range NN 
interaction that is dominated by pion exchange and constrained by chiral 
symmetry. The effects of the omitted Q-space, P{H^^f — H)P, might be 
expressed in some systematic expansion, with the coefficients of that ex- 
pansion directly determined from data, rather than from any knowledge of 
H acting outside of P. 

While we explored this second approach some years ago, some subtle 
issues arose, connected with properties of HO bases. This convinced us that 
the first step in our program should be solving and thoroughly understand- 
ing the effective interactions problem via the first approach, so that we 
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would have answers in hand to test the success of more phenomenological 
approaches. Thus we proceeded to follow the first approach using a real- 
istic NN potential, avlSj generating H^f^ numerically for the two- and 
three-body problems in a variety of HO SM spaces. Here I will use these 
rcsiilts to show that a properly defined short-range interaction provides an 
excellent representation of the effective interaction. This is an encourag- 
ing result, one that suggests a purely phenomenological treatment of the 
effective interaction might siiccecd. 

The key observation is that HOBET is an expansion around momenta 
A; ~ 1/6, and thus differs from EFT approaches that expand around k ^0. 
Consequently a HOBET P-space lacks both high-momentum components 
important to short-range NN interactions and long-wavelength components 
important to minimizing the kinetic energy. While our group has previously 
discussed some of the consequences of the combined infrared/ultraviolet 
problem in HOBET, here 1 identify another: a sharp energy dependence in 
H'^^-f that must be addressed before any simple representation of H'^^-f — H 
is possible. This, combined with a trick to remove operator mixing, leads to 
a simple and successful short-range expansion for H^-^^ — H.I conclude by 
noting how these results may set the stage for a successful determination 
of the HOBET H^'-ff directly from data. 

2. Review of the Bloch-Horowitz Equation 

The basis for the approach described here is the Bloch-Horowitz (BH) equa- 
tion, which generates a Hermitian, energy-dependent effective Hamiltonian, 
H^^-f , which operates in a finite Hilbert space from which high-energy HO 
Slater determinants are omitted: 



Here H is the bare Hamiltonian and E and are the exact eigenvalue 
and wave function (that is, the results of a full solution of the Schroedinger 
equation for H in P+Q). The BH equation must be solved self-consistently, 
as H'^f-I' depends on E. If this is done, the model-space calculation repro- 
duces the exact E, and the model-space wave function ^Pp is simply the 
restriction of \E' to P. If one takes for P a complete set of HO Slater deter- 
minants with HO energy < Apfkv, H^^f will be translational invariant. P 
is then defined by two parameters, Ap and the HO size parameter b. 

The BH equation was solved numerically for the avl8 potential using 



H^ff = H + H 




*p) = (l-g)|*). 



(1) 
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two numerical techniques. In work carried out in collaboration with C- 
L. Song,* calculations were done for the deuteron and '^He/^H by directly 
summing the effects of avl8 in the Q-space. Because this potential has a 
rather hard core, sums to 140 fkj were required to achieve ^ 1 keV accuracy 
in the dcutoron binding energy, and 70 to achieve ^ 10 kcV accuracies 
for '^Hc/''H. In work carried out with T. Luu,^ such cutoffs were removed 
by doing momentum-space integrations over all possible excitations. 

The results arc helpful not only to the goals discussed previously, but 
also in illustrating general properties of H'^^-f that may not be widely ap- 
preciated. For example. Table 1 gives the evolution of the P-space ^He avl8 
wave function as a function of increasing Ap for fixed &. \I/p evolves simply, 
with each increment of Ap adding new components to the wave function, 
while leaving previous components unchanged. One sees that the probabil- 
ity of residing in the model space grows slowly from its Ofuv value (31%) 
toward unity. 



Table 1. Evolution of the ^He avl8 HO wave function with Ap 



state 


amphtude 


Ohw 




4Itu> 


6fko 


8 Hlo 


exact 


(31.1%) 


(57.4%) 


(70.0%) 


(79.8%) 


(85.5%) 


(100%) 


1 0,1) 


0.5579 


0.5579 


0.5579 


0.5579 


0.5579 


0.5579 


|2,1> 


0.0000 


0.0463 


0.0461 


0.0462 


0.0462 


0.0463 


1 2,2} 


0.0000 


-0.4825 


-0.4824 


-0.4824 


-0.4824 


-0.4826 


1 2,3> 


0.0000 


0.0073 


0.0073 


0.0073 


0.0073 


0.0073 


|4,1> 


0.0000 


0.0000 


-0.0204 


-0.0204 


-0.0204 


-0.0205 


|4,2> 


0.0000 


0.0000 


0.1127 


0.1127 


0.1127 


0.1129 


|4,3> 


0.0000 


0.0000 


-0.0419 


-0.0420 


-0.0421 


-0.0423 



Effective operators are defined by 

and must be evaluated between wave functions ^p having the nontrivial 
normalization illustrated in Table 1 and determined by 

1 = (*|*) = (*p|i^//|f p). (3) 

The importance of this is illustrated in Fig. 1, where the elastic magnetic re- 
sponses for dciitcrium and ^Hc arc first evaluate with exact wave functions 
\E'p but bare operators, then re-evaluated with the appropriate effective 
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operators. Bare operators prove a disaster evcm at intermediate momentiim 
transfers of 2-3 f~^. By using the effective operator and effective wave func- 
tion appropriate to Ap, the correct result - the form factor is independent 
of the choice of Ap (or b) - is obtained, as it must in any correct apphcation 
of effective theory. 




q(fb 

Fig. 1. Deuteron and ^He elastic magnetic form factors evaluated for various P-spaces 

with bare operators (various dashed and dotted lines) and with the appropriate effective 
operators (all results converge to the solid lines). 

The choice of a HO basis excludes not only high- momentum components 

of wave functions connected with the hard core, but also low-momentum 
components connected with the proper asymptotic fall-off of the tail of 
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the wave functions. This combined ultraviolet /infrared problem was first 
explored by us in connection with the nonperturbative behavior of H'^^-f; 
the need to simultaneously correct for the missing long- and short-distance 
behavior of is the reason one cannot define a simple P that makes eval- 
uation of H''^-l^ converge rapidly. We also found a solution to this problem, a 
rewriting of the BH equation in which the relative kinetic energy operator is 
summed to all orders. This summation can be viewed as a transformation of 
a subset of the Slater determinants in P to incorporate the correct asymp- 
totic falloff. This soft physics, obtained from an infinite sum of high-energy 
HO states in Q, is the key to making H'^f^ perturbative. 

It turns out this physics is also central to the issue under discussion here, 
the existence of a simple representation for H'^^^ = H+H ^J^q^ QH, where 
H = T+V. The reorganized BH H'^f-^ is the sum of the three left- hand- side 
(LHS) terms in Eqs. (4) 

(a|T + rQ— Or|/3) — > {a\T\P) 

h/ — C^/i nonedge 
E E 

^ ^ E — TQ E — QT^^"^ nonedge ^ ^ 

^ ^ E — TQ E — QH^ E — QT^^'^ nonedge ^ ^ E — QH^ ^ ^ 

The first LHS term is the effective interaction for T, the relative kinetic 
energy. As QT acts as a ladder operator in the HO, E/E—QT is the identity 
except when it operates on an \a) with energy Apfko or (Ap — l)hiii. We will 
call these Slater determinants the edge states. For nonedge states, this new 
expression and the BH form given in Eq. (1) both reduce to the expressions 
on the right of Eqs. (4). 

Noting that the first LHS term in Eqs. (4) can be rewritten as 

E TOT E 

we see that the QT summation can be viewed as a transformation to a new 
basis for P, E/{E — QT)\a), that is orthogonal but not orthonormal. This 
edge-state basis builds in the proper asymptotic behavior governed by QT 
(free propagation) and the binding energy E. The transformation preserves 
translational invariance, as T is the relative kinetic energy operator. Viewed 
in the transformed basis, the appropriate effective interaction in given by 
the LHS terms in the square brackets in Eqs. (6) below. 

Alternatively, the results can be viewed as two equivalent expressions 
for the effective interaction between HO states, but with a different division 
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between "bare" and "rescattering" contributions 

E 



bare : ^ 



rescattering 



E-TQ 

E 



jj_TQT 



E 



E~TQ 



^^^^ 



E-QT 
E 



<^ H- — rr^QH (6) 



E-QT E-QH 



It is this new division that is critical. The expressions are identical for 
nonedge states. But for edge states, only the expression on the left isolates 
a quantity, VGV, that is short-range and nonperturbative. We will see that 
this is the term that can be represented by a simple, systematic expansion. 

Figure 2 shows the extended tail that is induced by E/E — QT acting 
on a HO state. Figure 3 is included to emphasize that there are important 
numerical differences between the two expressions in Eqs. (6). It compares 
calculations done for the deuteron using the two "bare" interactions: thus 
in both cases V enters only linearly between low-momentum states, and 
all multiple scattering of F in Q in ignored. Figure 3 gives the resulting 
deuteron binding energy as a function of b, for several values of Ap. For the 
standard form of the BH equation, a small model space overestimates the 
kinetic energy (too confining) and overestimates short-range contributions 
to V (too little freedom to create the needed wave- function "hole"). Making 
b larger to lower the kinetic energy exacerbates the short-range problem, 
and conversely. Thus the best 6 is a poor compromise that, even in a 10 /iw 
bare calculation, fails to bind the deuteron. But the new bare H on the LHS 
of Eqs. (6) sums QT to give the correct wave-function behavior at large r, 
independent of b. Then, for the choice b ~ 0.4-0.5 f, the short-range physics 
can be absorbed directly into the P space. The result is excellent Oth-order 
ground-state energy, with the residual effects of multiple scattering through 
QV being very small and perturbative.^ 

3. Harmonic-Oscillator-Based Effective Theory 

Now I turn to the question of whether (and how) the Q-space rescattering 
contribution to H'^^^ might be expressed through some systematic short- 
range expansion. There are two steps important in applying such an ex- 
pansion to HOBET. One has to do with the form of the short-range expan- 
sion. A contact-gradient (CG) expansion, constructed to include all possible 
LO (leading order), NLO (next-to- leading order) , NNLO (next-to-next-to- 
leading order), interactions is commonly used, 

afo(Ap, b)6ir) + a?,\o(Ap, &)( ^^^(r) + <5(r) V^) + 

<^N'NLoi^P, b)V'S{r)V' + a'j^'Zoi^p, b){V'S{r) + <5(r) V^). (7) 
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T 




r(fm) 

Fig. 2. A comparison of the \nl) and extended [E/E — QT)\nl) radial wave functions, 
for the edge state (n, I) = (6, 0) in a Ap = 10 deuteron calculation. Note that the 
normalization of the extended state has been adjusted to match that of \nl) at r=0, in 
order to show that the shapes differ only at large r. Thus the depletion of the extended 
state at small r is not apparent in this figure. 

Because HOBET is an expansion around a typical momentum scale 
^ l/b, rather than around fc = 0, it is helpful to redefine the derivatives 
appearing in the CG expansion Noting 

V"expifc-r |fe=o = 0,n= 1,2,...., (8) 
we demand by analogy in HOBET 

V"Vis(&) = 0,n- 1,2,... (9) 
This can be accomplished by redefining the operators O of Eqs. (7) by 

d e'-V20e'-V2 (10) 

The gradients in Eq. (7) then act on polynomials in r, a choice that 
removes all operator mixing. That is, if a^o is fixed in LO to the rt = 1 to 
n = 1 matrix element, where n is the nodal quantum number, it remains 
fixed in NLO, NNLO, etc. Furthermore, the expansion is in nodal quantum 
numbers, e.g., 

V2-(n-l) - (71- l)(n-2) (11) 
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0.4 0.6 0.8 1.0 1.2 

b(fm) 



Fig. 3. Dcuteron ground-state convergence in small P-spaces, omitting all effects due 
to the multiple scattering of V in Q. The standard BH formulation with P{T + V)P 
fails to bind the deuteron, even with Ap = 10. The reorganized BH equation, where QT 
has been summed to all orders but V still appears only linearly, reproduces the correct 
binding energy for Ap=6. 



SO that matrix elements become trivial to evaluate in any order. It can be 
shown that the leading order in n contribution agrees with the plane-wave 
result, and that operator coefficients are a generalization of standard Talmi 
integrals for nonlocal potentials, e.g., 

POC POO 

^NNLO- / e-^'^rlV{n,r^)rle-^'^rlrldndr, (12) 
Jo Jq 

The next step is to identify that quantity in the BH equation that should 
be identified with the CG expansion. This has to do with the two forms of 
the BH equation discussed previously. Consider the process of progressively 
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integrating out Q in favor of the CG expansion, beginning at A >> Ap and 
progressing to A = Ap. Using the projection operator 

A 

Qa= Qap=0 (13) 

a=Ap+l 

we can isolate the contributions, above some scale A, to the two BH rescat- 
tering terms of Eqs. (6) 



E 

Aot(A) 



E-TQ 



^ .(14) 



E-QT 



The goal of a CG expansion might be successful reproduction of the matrix 
elements of A (A) - the Q-space rescattering contributions for the standard 
form of the BH equation - as A — > Ap. This would allow us to replace all Q- 
space rescattering by a systematic short-range expansion, opening the door 
to a purely phenomenological determination of H'^f^ for the SM. The test 
case will be an Shw P-space calculation for the deuteron (Ap = 8, 6=1.7 
f). The running of the 15 independent ^Si matrix elements of A(A) are 
plottted in Fig. 4a. Five of these are distinguished because they involve an 
edge-state bra or ket (or both). The evolution of these contributions with 
A is seen to be somewhat less regular than that of nonedge-state matrix 
elements. The results for A = Ap show that rescattering is responsible for 
typically 12 MeV of binding energy. 

The CG fit to the results in Fig. 4a were done in LO, NLO, and NNLO as 
a function of A, using the standard form of the BH equation. The coefficients 
are fit to the lowest-energy matrix elements. Thus in LO aLo{^) is fixed 
by the Is — Is matrix element, leaving 14 unconstrained matrix elements; 
the NNLO fit (Is — Is, Is — 2s, Is — 3s, and 2s — 2s) leaves 11 matrix 
elements unconstrained. This is easily done, because the operators do not 
mix; e.g., among these four, only the Is — 3s matrix element is influenced 
by a^NNLO- result is a set of coefficient that run as a function of A 
in the usual way, with qlo small and dominant for large A, and with the 
NLO and NNLO terms turning on as the scale is dropped. Figs. 4b-d show 
the residuals the differences between the predicted and calculated matrix 
elements. For non-edge-state matrix elements the scale at which typical 
residuals in A are significant, say greater than 100 keV (above ~ 1%), is 
brought down successively, e.g., from lOOhuj, to ^ 55Hlu (LO), to ^ 25huj 
(NLO), and finally to ~ Aphiu (NNLO). But matrix elements involving 
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edge states break this pattern: the improvement is not significant, with 
noticeable deviations remaining at ~ lOOSw even at NNLO. 




Fig. 4. In a) rcscattcring contributions to H'^ff due to excitations in Q above A are 
given for the standard form of the BH equation. These results are for the 15 matrix 
elements that arise in an 8?iai calculation for the deuteron, with 6=1.7 f. In b)-d) the 
residuals of LO, NLO, and NNLO fits are shown (see text). Matrix elements with bra or 
ket (dashed) or both (dot-dashed) edge states are seen not to improve systematically. 



This failure could be anticipated: because QT strongly couples nearest 
shells across the P — Q botmdary, H ^^^j^ QH contains long-range physics. 
The candidate short-range interaction is V ^\jj QV , not QH: this 
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is the reason QT should be first summed, putting the BH equation in a form 
- the LHS of Eqs. (6) - that isolates this quantity. This reorganization 
affects edge-state matrix elements only, those with the large residuals in 
Figs. 4b-d. 

To Msc the reorganized BH equation, the QT sums appearing in Eqs. 
(4) must be completed. There are several procedures for doing this, but 
one convenient method exploits the raising/lowering properties of T. The 
result is a series of continued fractions gi{2E/tiuj, {ai}, {A}), where = 
(2n+2i+Z-l/2)/2 and/3i = ^{n + i){n + i + I + l/2)/2. For any operator 
6 (e.g., V, V^^QV, etc.) 

^'''^'^Y^^Y^^''^^ = ,g,(n',r)ff,(n,Z)(n'+j l\0\n + i I) (15) 

It follows that the coefficients of the CG expansion for a HO basis must be 
redefined for edge states, with a state- and E-dependent renormalization 

aLO^a'i^Q{E/hio;n',l',n,l) = aLO ^ Ojin' ,l')gi{n,l) 

[ T{n' + j + l/2)r{n + i + 1/2) 
^[ r(n' + l/2)r(n+l/2) 

This renormalization, which introduces no new parameters, can be evalu- 
ated in a similar way for heavier systems: T remains a raising operator. 

Figs. 5 shows the results: the difficulties encountered for A (A) do not 
arise for Aqt(A). The edge-state matrix elements are now well behaved, 
and the improvement from LO to NLO to NNLO is systematic in all cases. 
When A — > Ap, the CG potential continues to reproduce H^^^ for the avl8 
potential remarkably well, with an the average error in ^Si NNLO matrix 
elements of about 100 keV (or 1% accuracy). Other channels we explored 
behaved even better: the average error for the 15 ^Sq matrix elements is 
about 10 keV (or 0.1% accuracy). Because all matrix elements of H^^^ are 
reproduced well, the CG potential preserves spectral properties, not simply 
properties of the lowest energy states within P. The NNLO calculation in 
the ^Si-^Di channel yields a deuteron ground-state energy of -2.21 MeV. 

Several points can be made: 

• The net effect of the QT summation is to weaken the CG potential for HO 
edge states: the resulting, more extended state has a reduced probability at 
small r. Consequently the effects of QV are weaker than in states immune 
from the effects of QT. 

• The very strong QT coupling of the P and Q spaces is clearly problematic 



J./ ^ 



(n'- l)!(n- 1)! 
(n' + j - l)\{n + i - 1)! 



(16) 
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Fig. 5. As in Fig. 4, but with the edge states treated aceording to the QT-summed 
reorganization of the BH equation, as described in the text. 



for an ET; small changes in energy denominators alter the induced inter- 
actions. Thus it is quite reasonable that removal of this coupling leads to a 
strong energy dependence in the effective interaction between HO states. I 
believe that proper treatment of this energy dependence will be crucial to 
a correct description of the bound-state spectrum in the HO SM. 
• This process can also be viewed as a transformation to a new, orthogonal 
(but not normalized) basis for P in which la) e-^qt I'^)' This yields basis 
states with the proper asymptotic behavior for each channel. A CG expan- 
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sion with fixed coefficients can be used between these states, following the 
reorganized BH equation of Eqs. (4). 

• While our calculations have been limited to the deuteron, this same phe- 
nomena must arise in heavier systems treated in HO bases — QT remains 
the ladder operator. As the issue is extended states that minimize the ki- 
netic energy, it is clear that the relevant parameter must be the Jacobi 
coordinate associated with the lowest breakup channel. This could be an 
issue for treatments of ^ based on the Lee-Suzuki transformation, which 
transforms the interaction into an energy-independent one . In approaches 
like the no-core shell model, the Lee-Suzuki transformation is generally 
not evaluated exactly, but instead only at the two-body level. If such an 
approach were applied, for example, to ^Li. a system weakly bound (1.475 
MeV) to breaking up as a-l-d, it is not obvious that a two-body Lee-Suzuki 
transformation would treat the relevant Jacobi coordinate responsible for 
the dominant energy dependence. This should be explored. 

• I believe the conclusions about the CG expansion will apply to other ef- 
fective interactions. For example, V-low-k,^^ a soft potential obtained by 
integrating out high-momentum states, is derived in a plane-wave basis, 
where T is diagonal. Thus it should be analogous to our CG interaction, 
requiring a similar renormalization when embedded into a HO SM space. 
It would be interesting to test this conclusion. 

4. Summary 

These results show that the effective interaction in the HO SM must have 
a very sharp dependence on the binding energy, defined as the energy of 
the bound state relative to the first open channel. This is typically to 10 
MeV for the bound states of most nuclei (and 2.22 MeV for the deuteron 
ground state explored here). Once this energy dependence is identified, the 
set of effective interaction matrix elements can be represented quite well 
by a CG expansion, and the results for successive LO, NLO, and NNLO 
calculations improve systematically. 

This result suggests that the explicit energy dependence of the BH equa- 
tion is almost entirely due to QT - though this inference, based on the 
behavior of matrix elements between states with different number of HO 
quanta, must be tested in a case where multiple bound states exist. 

We also presented a simple redefinition of the gradients associated with 
CG expansions, viewing the expansion as one around a momentum scale 
1/6. This definition removes operator mixing, making NNLO and higher- 
order fits very simple. The expansion then becomes one in nodal quantum 
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numbers, with the coefRcients of the expansion related to Talmi integrals, 
generalized for nonlocal interactions. 

While our exploration here has been based on "data" obtained from 
an exact BH calculation of the effective interaction for the a.vl8 potential, 
this raises the question, is such a potential necessary to the SM? That is, 
now that the success of an NNLO description of Q-space contributions is 
established, could one start with PHP and determine the coefficients for 
such a potential directly from data, without knowledge of matrix elements of 
H outside of P? I believe the answer is yes, even in cases (like the deuteron) 
when insufficient information is available from bound states. It turns out 
that the techniques described here can be extended into the continuum, 
so that observables like phase shifts could be combined with bound-state 
information to determine the coefficients of such an expansion. An effort of 
this sort is in progress. 

I thank M. Savage for helpful discussions, and T. Luu and C.-L. Song for 
enjoyable collaborations. This work was supported by the U.S. Department 
of Energy Division of Nuclear Physics and by DOE SciDAC grant DE- 
FG02-00ER-41132. 
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